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TECHNICAL MEMORANDUM NO. 1112 

THE EFFECT OF TURBULENCE ON THE FLAME VELOCITY IN GAS MIXTURES! 

By Gerhard Damkohler 


The present report deals with the effect of turbulence on the propa- 
gation of the flame. Being based upon experiments with laminar as well 
as turbulent Bunsen flames, both the physico-chemical and the hydro- 
dynamical aspects of the problem are analyzed. 

A number of new deductions, interesting from the point of view of 
engine combustion and other very rapidly changing flame reactions, are 
made. 


SUMMARY 


1. ; By use of the Bunsen burner method on various propane-oxygen mix- 
tures, a number of flame velocity measurements were made in the range of 
Reynolds numbers Re= 612 to 17300 (referred to inside diameter of burner 
tube), so that the appearance of the turbulent Bunsen flame could be com- 
pared with that of the laminar. Tho inside diameters d of the employed 
tubes I, II, and III were 1 . 385 * 2.18, and 2.718 millimeter^ respectively. 

2. The laminar Bunsen flame burned completely without noise. It 
showed a sharply defined bell— shape flame surface, the base of which ex- 
tended a little over the inside cross section of the burner tube. The tip 
of the cone was rounded off with a curvature radius of the same order of 
magnitude as the luminous thickness Or, of the flame surface itself; the 
thickness of the luminous zone 5-^ was measured from time to time just a 

little below the cone tip and amounted to 0.10 -*0.02 millimeter in all 
the laminar tests. 

3. The turbulent flame was recognized acoustically by loud whistling. 
The zone of the flame, which in its base, also extended over the internal 
cross section of the burner tube, was no longer sharply defined, but 


i"Der Einfluss der Turbulenz auf die Flammengeschwindigkeit in 
Gasgemischen . " Zeitschrift fur Elektrochemie und nngewandte 
Physikalische Chemie. 'vol. 46, no. 11, Nov. 1940, pp. 601 — 626 . 
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washed out. The zones were no longer hell-shape except on the narrow- 
est tube I; they were already rather holt-shape on the wider tube II, 
and on tube III its upper part was considerably enlarged with respect 
to its base. 

4. From the boundary surfaces of the turbulent flame zone, even 
though blurred, the existing maximum and minimum flame velocities, 

Wf ma;x; and w f man , can be approximately computed. 

With exception of two tests the minimum flame speeds w .■ m |. n 
agreed consistently with the corresponding flame speed w^q for laminar 
flow. The two unsuccessful tests can, however, be interpreted likewise. 


II with Re - 2300 to 5000) approximately proportionally to v Re o.t large 
Reynolds numbers (tubes II and III with Re = 5000 to 18000) approxi- 
mately linearly with Re. 

5. The theory of laminar, stationary flame zone is briefly disr- 
cussed. The terms "preparation zone", "reaction zone", "ignition tem- 
perature", "rate of chemical reaction", and "flame speed, "are explained, 
and the mathematical relations existing between these quant 1 ties de- 
duced. Numerical values for the quantities are given on the basis of 
the experimental material. 

6. The mean chemical decomposition velocity in the reaction zone 
amounts to about 0.2 mol of fuel per cubic centimeter times second? 
in the tests, so that every fuel molecule receives approximately 10 b 
impulses before being completely burned. 

7. Comparison of the computed thicknesses of the preparation and 
reaction zones, f y and. S.p, respectively, with the experimentally 
established luminous zone thickness 6^, gives the fairly arbitrary 

alternative: the ignition temperature « z ~ 500 C in the flame either 

agrees approximately with the autoignition temperatures determined for 
the particular mixture by the method of adiabatic compression or with 
the Dixon burner arrangement; in which case 5^ ~ 6y + 5^; or else with 

t 200° C the ignition temperature in the flame lies considerably 
lower than the cited autoignition temperatures; in which case 5^ ~ . 

The second possibility appears more likely at the present time. Inves- — 
tigatien of the intensity distribution of the individual emission spec— 
trums in 5-^ would alao afford insight into the time sequence of it3 
occurrence and thus contribute considerably to the elucidation of the 
chemical reaction mechanism in the flame. 


At small Reynolds numbers the ratio 


I 


f max 
w f min 


increases (tube I and 
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8. The propagation of a laminar free jet is numerically estimated 
after an article ty H. Schlichting and found to be compeletely negligi- 
ble for the interpretation of the laminar Bunsen cone. However, the 
experimentally observed laminar flame cones do not completely agree in 
form with cones computable on the assumption of a parabolic velocity 
profile and equal flame speed at every point of the flame cone. The 
discrepancies can be explained. Directly above the burner tube rim a 
horizontal, outwardly directed flow must exist, which can be numeri- 
cally appraised, and which effectively prevents air from the outside 
being diffused into the cone of the flame. 

9. Prandtl r s characterization of the turbulence of a flow by two 

numerical quantities is briefly outlined: by the mixing length Z, 

which indicates the approximate magnitude of the turbulence bodies, and 
by the exchange quantity e, which, as a type of diffusion coefficient, 
permits a quantitative statement concerning the amplified transport 
processes in turbulent flow. With these two quantities the turbulence 
in the free jet issuing from a tube orifice is numerically estimated, 

it itself being divided into a nuclear zone with the flow and turbu- 
lence properties as in a simple pipe (explored by Wikuradse) and a bor- 
der zone with particular characteristics (investigated by W. Tollmien) . 

10. The theoretically conceivable influences of the turbulence on 

the propagation of flames are discussed. Two types of turbulence with 
totally different effects must be distinguished: "coarse— body turbu- 

lence with a great mixing length compared to the thickness of the lami- 
nar combustion zone, with roughened and ultimately scattered flame 
front, in which, as a result, the surface of the flame is enlarged, and 
the flame speed referred to the flow cross section rises beyond the 
laminar value Wj*j_. On the extremely "fine— body" turbulence with very 
short mixing length compared to laminar flame zone thickness, there is 
no roughening of flame surface. On the other hand, all transport pro- 
cesses in it are greater; hence also the rise in turbulent flame speed 
w ft beyond the laminar w.pp , which is, according to the formula, 



The values of T?ft computed with it (averaged over the total flow 
cross section) are only about 30 percent higher than the experimen- 
tally established maximum flame speeds in the turbulent flame. Hence 
in the tests, besides the rough— * body turbulence, the fine— body turbu 
lence also must have played a part; especially since the greatest mix- 
ing lengths in the three burner tubes amounted to 0.11, 0 . 17 , and 
0.21 millimeter respectively, while the luminous zone thickness was 
fixed at 8^= 0.1 millimeter. 
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For the combustion in the Otto-cycle engine the coaree-body turbu- 
lence is likely to be decisive, although its effect is as yet unpredict- 
able; because turbulence bodies of the order of Z =0.1 millimeter can 
be produced only for gas issuing from very narrow flow cross sections, 
and such are always to be avoided for reasons of the high flow resist- 
ances involved. 

11. In turbulent interchange or mixing processes it is always nec- 
essary to distinguish between transport lengths that are greater or 
smaller than the mixing length Z. In the latter case the Prandtl inter- 
change quantity is effective only with a fractional amount of its macro- 
scopic value e. 


INTRODUCTION 


The earlier reports on flame speed (ignition velocity) in gas mix-, 
tures deal exclusively with the propagation of flames in static or lam- 
inar gas flow (in the bomb, in the pipe closed, at both ends, in the soEp 
bubble, or in the Bunsen burner, etc). The reports further define the 
physical and chemical influences (pressure, initial temperature, type 
of container walls, chemical composition of initial gas) and attempt 
to explain the phenomenon of flame propagation theoretically and from 
it deduce physical and kinetic reaction data. 

But not one of these reports contains any information about the 
propagation of flames in turbulent gas flow, although this case is of 
considerable importance for the technique. In the Ottc-cycle engine, 
for instance, flame speeds of from 10 to 50 meters per second occur, 
which are due solely to turbulence effect (measured relative to cylinder 
head), that is, from 10 to 100 times the amount of the value measured 
on the same gas mixture without turbulence. However, it still remains 
to be explained whether these high combustion speeds in the engine are 
simulated by gas masses that flow faster than they burn or whether real- 
ly higher flame speeds are involved, the cause of which might be found 
in the turbulent motion, but likely also in dynamic gas effects. 

The little knowledge of the influence of pure turbulence on the 
flame speed may be due to the following causes: in the first place 

it is difficult to identify the turbulence of a flow quantitatively. 

So far this has been successfully accomplished only for the developed 
turbulence in stationary flow. (Cf. the fundamental works by L. 

Prandtl and his school as well as the more recent British and American 
wind-tunnel investigations on the subject. ) On top of that these re- 
ports are scarcely known to engineers and physico-chemists interested 
in the problem of flame speed. 
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The subsequent report is a first attempt to measure exactly the in- 
fluence of turbulence on f lame propagation and to interpret it theoreti- 
cally. To this end comparable flame speed measurements were made in 
laminar and fully developed turbulent flow on various uniformly flowing 
gas mixtures. It is true that the latter case corresponds only in part, 
to the conditions encountered in the Otto engine, but it also makes it 
possible to give numerical data on the state of turbulence of the flow 
which are indispensible for a deeper insight into the processes involved. 

The theoretical treatment of the experimental data begins with the 
development of the general theory of the stationary flame zone in lami- 
nar flow; then f ollow the numerical data on preparation zone, reaction 
zone, ignition temperature, and chemical reaction speed in the cone of 
the laminar Bunsen burner, the enlargement of a laminar gas jet issuing 
into the open and the form of Buncen cone burning in it, quantitative 
appraisal of the conditions of turbulence in the free jet in the light > 
of the Prandtl school, discussion of the theoretically conceivable tur- 
bulence effects on the flame speed, and lastly, comparison with the 
test data of the present study. 


TEST METHOD AND TEST DATA 
1. Brine ip les 


The measurement of the effect of turbulence on flame propagation 
in a flowing gas mixture is predicated upon the numerical data on the 
turbulence . 

The only method of flame speed measurement in which definite and 
well- known flow conditions exist and in which wall effects are also ex- 
cluded as far as possible, is the Bunsen burner method. Here a flame 
cone stands stationary in an open jet that issues from a burner tube 
and in its core still manifests the same turbulence as the flow in the 
tube. Very carefully executed measurements of turbulence in pipe flow 
are available (reference 1); while the turbulence in the boundary zone . 
of the free jet, as investigated mathematically by V/. Tollmien (refer- 
ence 2), showed good agreement with Gottingen wind-tunnel measurements. 
Through these studies the turbulence in the free jet is known. 

Admittedly, the simple burner tube has the drawback of a locally 
dependent flow velocity and turbulence, hence necessitates integration 
over the entire flow cross section. But then the turbulence is known 
and need not be measured first. 
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2. Experimental Setup 

The setup, is diagrammatic ally shown in figure 1. The gases from 
storage tanks Aj. and A 2 pass, at first, separately across the bubble 
valves Bi and B 2 (with controlled immersion depth) filled with water or 
mercury, depending upon the flow velocities, into the calibrated flow 
meters sti and st 2 . From there they pass across the pressure- equal- 
izing flasks Dj. and D 2 and the reducing valves Vj. and V 2 into the mixing 
piece T, topped by the actual burner tube B. Pressure fluctuations wore 
avoided by the control of valves Vj. and V 2 in conjunction with flasks 
Di and D 2 ; hence the flame burned perfectly still. A glass tube R (80 
mm in diameter, 300 iron long) slipped over the burner provided protection 
against draft. The dimensions of the brass burner tubes are given in 
table 1. 


TABLE 1 

DIMENSIONS OF BURNER TUBES 


Tube 

Length 

Diameter (mm) 

Wall thickness L 

at mouth d 

Outside i 

1 

i 

i Inside 

! 

1 

I 

500 

j 

2.0 

! 

1.385 

0.30 

361 

II 

500 

4.0 

2.18 

Ml 

229 

III 

1000 

5.0 

i 

2.718 

L 

' .60 

368 


Tubes: I find II were tapered at the upper end over a 20-millimeter 
length so that the outside diameters and hence the wall thicknesses de- 
creased to the figures cited in the table. The relative entrance 
lengths g- exceed the value of 150 and 50 necessary to produce devel- 
oped f low several times over in laminar flow (reference 3) as well as 
in turbulent flow (reference 1). These lengths were intentionally 
chosen greater in order to ensure complete mixing of the gases, which 
was further aided by the T piece. In the preliminary tests the flame 
cones necessary for the calculation of the flame speed were measured 
direct with the small telescope of a cathetometer, in the principal 
tests by photographic recording and subsequent measuring of the positive 
prints. A scale near the burner orifice (fig.: 1) was also photographed 
(1 division = l/2 mm). After considerable trials the Afga— Kontraet 
plate 6.5 X 9 square centimeters proved the most suitable. 
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3 . Experimental Gases 

The gases consisted of various blends of technical oxygen and 
technical propane, which is essentially a propane-butane mixture, con- 
taining about 6 percent of unsaturated hydrocarbons (olefins). It 
might be suspected that more rapid vaporization of the easier boiling 
propane in the original steel flask would gradually enrich the rest of 
the butane; but this was not the case as proved by several molar weight 
determinations!, which, while manifesting fluctuations, indicated no 
systematic time process. It is likely that the inferior diffusion coef— 
ficient of the fluid propane— butane mixture provided for practically 
total vaporization on the surface of the liquid. The mean mole weight 
of propane in the tests was M = 46.1 £ 0.5 correspoinding to a propane 
content of 85 — 5 percent, if the rest is regarded as butane. 

The Reynolds number in the burner tube was computed by moans of 
the viscosity of the gas mixtures with the apparatus shown in figure 2. 
It consisted largely of a glass tube spiral of about 1-millimeter inside 
diameter. The entrance length . AB was about 500 millimeters, the actual 
test length BC, 1000 millimeters. The flow is laminar at Reynolds 
numbers ranging between 100 and 500 , hence the following relation for 
pressure drop between B and C; 


whore 


Ap = k r) G 


r, dynamic viscosity 
G flow volume 
k • constant 

The latter was empirically determined by calibration measurements with 
pure oxygen at 47xl3xl0 4 ,with t) 0 ^ = 2.05^ xlCr 4 grams per centi- 
meter -1 second "* 1 . 

The recorded viscosities are reproduced in table 2 and figure 3> 
they were all made at l 8 + 1° C with a propane of density 

1.92 x 10“3 grams per cubic centimeter (M = 45.9)* 

4. Experimental Results 

The subsequent data refer almost exclusively to the principal 
tests, because in the preliminary tests in which the flame cones were 


1 With the Kahle type gas sweep obtained from Linde's Eismuschinen 
(Hollriegelskreuth, near Munich) . 
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measured visually with eathe tome ter, the extent of dispersion of the 
observed flame velocities amounted bo 20 percent, as against less than 
10 percent on the averaged main test curves (fig. 5)* Particularly 
uncertain in the preliminary tests was the gaging of the blurred turbu- 
lent flame zones, since the field of vision of the small telescope cov- 
ered only part of the flame; whereas the positive print of the pho bo- 
graphs in the principal tests always presented the entire flame end made 
it much easier to identify individual zones. These zones were then out- 
lined by pencil and measured with the telescope of the cathetometer. 

The photographs of the individual flames are reproduced in figure 4 and 
tabulated in tab le 3. The flame phenomena obtained were distinctly 
different, depending upon whether the flow in the burner tube was lami- 
nar or turbulent. 

The laminar flow (Re« 2300) burnt absolutely without noise. It 
revealed a sharply defined bell— shape f lame surface, the base of which 
extended over the internal cross section of the burner tube. The tip 
of the cone was rounded off with a curvature radius of the same order 
of magnitude as the thickness of tho luminous flame surface. This 
luminous cone thickness was measured from time to time j.'uet a little be- 
low the tip of the cone and ranged between 0. lO^O.OP milliia-jter through- 
out all laminar tests. 

The turbulent flame (Re » 2300) was accompanied by a loud 

whistle. The zone of flame, which in its base also extended over the 
internal cross section of the burner tube, was no longer sharply de- 
fined but blurred. Only the narrowest tube I disclosed a bell-shape 
form; on tube II it became bolt-shape and on tube III' the upper part 
was already c ons i dor ab ly wider than at the base. 

For the quantitative interpretation of the test data the flame 
velocity wf (cm/sec) averaged over the flame cone was calculated; it in- 
dicates the speed at which the flame burns against the stationary gas 
flow. According to Gouy the following relation (reference 4) is 
applicable : 

w f = G/F (1) 

where 

G amount of gas supplied to "burner, cubic centimeters' per second 
F surface of flame cone, square centimeters 

With F computed as the surface of a body of rotation by Passus 1 . ^ 

theorem F = 2jtR t S from the length S of the enveloping line of the 
combustion cone outlined by pencil end the distance E * of its center 
of gravity from the center line, the outer or inner boundary surface 
of the flame zone can be used as a basis. The former served in the 
preliminary tests because it was easier to measure by the visual 
observation of the flame, while in the principal tests one flame 
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velocity each was determined for both luminous zone boundaries, the in- 
tention being to evaluate laminar and turbulent tests as nearly alike 
as possible. Closer examination of figure 4, especially the comparison 
of photographs 49 and 50 with 51 and 53 reveals that the turbulent flame 
zone issues from the laminar in such a manner that the boundary surfaces 
of the laminar luminous zone shift farther apart and at the same time be- 
come hazy. But in the turbulent flame these fairly distant and, in 
addition, blurred boundary surfaces are the only definite areas that can 
be measured and both of them must therefore be considered if a somewhat 
quantitative picture of the total flame phenomenon is to be secured. 


Both boundaries have, moreover, a physical significance also; the 
inner is the locus of the speediest combustion, the outer, of the slow- 
est combustion. The slowest flame velocities of the turbulent flame 
for all three tubes agreed with the flame velocity for laminar flow 
(identical mixture presumed) according to table 3 and figure 5* This 
result is very obvious; for a flame velocity slower than for laminar 
flow can be visualized only when the initial gas, before its complete 
combustion is di Luted by incombustible gas additions (for instance, 
products of combustion), which would be conceivable in the upper parts 
of a flame for strongly turbulent mixing motion. This condition prob- 
ably takes part in tests 52 and 45 which alone fell out of line. 
Possibly the "marginal turbulence" of the free jet exerts here such a 
delaying effect. In test 52 wf m i n is smaller than the corresponding 
laminar flame velocity (cf . fig. 5) and test 45 disclosed an irregu- 

, . j . , , , . v f max , _ „ . > 

larity m the ratio (cf. fig. 6). 

w f min _ 

„ . ,, , . w f max . „ , , 

So, since the ratio is a measure for the increase 

w f min 


in flame velocity due to turbulence, it was plotted in figure 6(a) and 


6(b) against the Reynolds number Re = — of the tube flow and against 
D 

-/"Re. The graphs show: for 2300< Re< 5000 (1.385 and 2.l8 mm tube 

diameter) , - increases approximately propcrtionately to ^/Re,' 


w f min 

for 5000 < Re < 18000 (tube diameters 2.18 and 2.78 mm) approximately 
linearly with Re. For the present, these experimentally obtained re- 
lations should be regarded as approximate formulas, because of the mar ad 
scattering of the test points. 


Lastly, it may be mentioned that when using a burner tube of only 
500-millimeter length and 2. 718— millimeter diameter, hence with a 
relative entrance length L/d = 184, the outer boundary of the turbu- 
lent flame zone revealed from '‘ime to time single, plainly distinguish- 
able gas bodies, which burned by themselves. But this phenomenon was not 
repeated on the twice-as-ldng tube III. Therefore the entrance length 
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of the 500~miilimeter tube was apparently insufficient for complete mix- 
ture of the original gases. 


III. THEORY - DISCUSSION OF TEST DATA 
1. Theory of Stationary Flame in Laminar Flow 


The flame front is divided into a preparatory zone (in the litera- 
ture usually termed "preheating zone") (refererce p) and a reaction zone. 
At the boundary surface of these two zones the ignition temperature {t z 
prevails. This demarcation is also useful as approximation for the mod- 
ern conception of the combustion processes; for these always have, as 
chain reactions, a certain induction period, during which the mixture, 
by forming or diffusion of active intermediate products as well as by 
heat input, must first be prepared for ignition. Here the entire chemi- 
cal decomposition remains at first quite email, usually less than 1 per- 
cent. At the place of the ignition temperature the concentrations 

of the active chain carriers and the reaction velocities of the individ- 
ual elementary reactions have by then become so great that the total de- 
composition velocity suddenly speeds up: here the reaction zone starts. 


Only in one respect do the modern conceptions show an essential 
difference from the earlier ones: the ignition temperature may not be 
regarded as a characteristic physical constant of the initial mixture, 
as pointed out by Jost and Von Muffling (reference 6). Neither can the 
ignition temperature i> z occurring in the flame on the border between 
preparatory zone and reaction zone be expected to agree with the igni- 
tion temperatures that are obtained for the same initial mixture by the 
Nernst method of adiabatic compression or Dixon’s burner tube. It is 
likely that the ignition temperature v z in the flame is lower, because 
ignition takes place under more favorable conditions than in the other 
two cases, in which the intermediate products necessary for ignition 
must be formed pure ly thermally. In the preparatory cone of the flame, 
on the other hand, these substances are supplied by diffusion from the 
reaction zone. Thus the new thermal formation of active substances can 
be smaller or be absent altogether. That the diffusion flow of these 
active intermediate products against the actual gas flow is substantial, 
has been demonstrated by Jost and Von Muffling. 


The conditions in the flame front in stationary flow are defined by 
the following equation system (reference 7 ) • . Conservation of momentum: 


P (w v) w = 


yp+T)[A w + - grad div w] 


( 2 )- 
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Conservation of energy: 

0 = QU - div (cp pwd div (-X gradd) (3) 

Conservation of atom type : 

0 = -Vj U - div (n 1 w). - div (-Bj grad n^) (4) 


where : 

w flow velocity 

P density 

p pressure 

£ temperature 

c p specific heat of mixture 

■n dynamic viscosity 

X thermal conductivity 

Q heat of reaction 

U chemical reaction speed (mol/cm 3 sec) 

n. number of mols of the particle type j per unit of volume 
J 

u j stoichiometric conversion factor 

D * diffusion coefficient of the particle type j 

The solution requires the integration in the preparatory zone with 
U = 0 and in the reaction zone with U ^ 0. 

Assuming the same total pressure at all points (which holds good 
accurately enough in the slow combustion of normal flames, but not in 
detonation processes, of course) and limited to one— dimensional flow, 
equation (2) cancels and w is simply defined by the equation of con— • 
tinuity 


P w = constant = P 
— a -a 


(5) 


subscript a denoting the initial state. 
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lastly, there is the equation defining the laminar flame velocity 

Yfl = Istal (6) 

Integration of equation (3) over the preparatory zone (U = 0) be- 
tween the initial state .at x = — <» and the point of ignition (subscript 
z) at x = 0, where X— - = 0 for x = — 00 , gives 

dx 

? p p a K a ( * z ~ 3a ) « X z )z (7) 


with Cp denoting the mean specific heat of flowing gas between -9 a 
and d z . Assuming the thermal conductivity in the preparatory zone as 
being constant, or computing with a suitable average value X = k, the 

integration of equation (3) gives the temperature variation in the pre- 
paratory zone : c p P a w a x 

fl ~*a = (*z- «a ) e X (8) 


(also given in reference 4, p. 20) . From this a thickness of the pre- 
paratory zone 5v = — x can be estimated when assuming for this x, 

for example: * — ~ - a = — 1~. 

■9 z —'9 a 100 

It then yields 


S p a-a (“Sv) 


lg 


lo 


(9) 


or, with allowance for equation (6) : 

5X 


&v~ 


c p p a w fl 


( 10 ) 


The thickness of the preparatory zone is thus reduced to the deter- 
mination of the laminar flame velocity wpp. Unfortunately it has not 
been possible up to now to check the thus computed values of 6 V by 
actual experiment. But they will be useful for the understanding of the 
turbulence effect on the flame velocity. 

An integration of equation (4) (there is one equation each ‘for 
each type of molecule) over the preparatory zone would result in rela- 
tions similar to equation (7), except for the concentrations nj sub- 
stituting for temperature -9 and the diffusion coefficients Dj re- 
placing the thermal conductivity X. Several additional assumptions 
would even afford equations corresponding to equation (8), and which 
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would describe the penetration of the active intermediate products fro:., 
the reaction zone into the preparatory zone, so far as a new thermal 
formation cf these intermediate products in the preparatory zone itself 
can be excluded. 

•Important xor xne suosequent study is the integration cf that - Ra- 
tion (4) which refers to the fuel particle itself (subscript b instead 
of j), the integrations being directly extended to cover the entire 
flame zone, or preparatory zone + reaction zone. Now , if the stoichio- 
metric conversion factor for the fuel particle is standardized to 

= 1, that is, if U indicates the number of molsof fuel per volume per 
unit time converted at a point of the reaction zone (with U - 0 in the 
preparatory zone as before), and it is borne in mind that at both limits 
x - ~o(x~ — Sy, respectively), and x = Sp (end of reaction zone) the ex— 

on- n 

pressions 3> 0 — - disappear, and further that n^ = at start and 

nh ~ 0 at finish, 

U X 8 R (11) 


( 12 ) 

U represents the mean chemical reaction speed (mol/cn 3 sec fuel) in the 
reaction zone. Hence U is directly United with the thickness of the 
reaction zone S R and the laminar ignition velocity w^. Naturally, 

equation (12) could also have been posted immediately for continuity 
reasons, but this method affords a better picture of the physical view 
points that play a part in the solution of the complicated system of equa- 
tions (3) and (4) by approximation. In the preparatory zone it was the 
temperature rise 1 ; in the reaction zone it was the disappearance of the 
fuel with a mean cnemical conversion speed U. 

The relation ensuring the continuous connection between reaction 
zone and preparatory zone is obtained as follower putting, similar 
to A. Eucken (reference 8) 



n-ba H a = 


or with equation (6) : 


J 


r 


&R 

U dx = 


r 


OR 

U dx = 


L 


o 

n ba vfi = U 5 r 


!For which reason the earlier reports always spoke only of the pre- 
heating zone, whereas the term "preparatory zone" fits the nature of the 
phenomenon much better. 


14 


MCA TM Nc. 1112 


f signifying, for the time being, an open numerical factor, but which 
cannot be overly far from unity, and ■ 3 e deno'ting the terminal combustion 
temperature. This becomes apparent on the actual temperature curve in 
the flame zone, which is plotted to scale in figure 7 (a) and 7 (t) on the 
basis of subsequently computed numerical values. From ( 6 ), (7) and (13 ) 
the thickness of the reaction zone follows 


c p Pa w fl(dz —'3a) 


(14) 


Elimination of from (12) and (14) leaves essentially a relation 

between flame velocity w^^, thermal conductivity X z and mean chemical 
reaction speed U: 


fX z U (d e - -3 Z ) 

c p r *ba ~ ‘3a) 


(15) 


Naturally, equation (15) can also be solved with respect to U; hence IT 
in the laminar flame can be obtained from experimentally accessible quan- 
tities. This will be done later, at which time it will also be consid- 
ered whether U in the normal flame is actually so small that — from the 
pure kinetic reaction point of view — still substantially higher U 
va.lues can be visualized, which should appear then, for example, in deto- 
nation processes. 


Formula (15) is essentially agreeable with the relations of Jouget, 
Nusselt and Daniel, but rather than assume a special law for the chemical 
reaction speed U, the factor f is introduced. 


2. Numerical V alues for Preparatory Zone, Reaction Zone, Ignition 
Temperature and Chemical Reaction Speed in the Cone 
of the Laminar Bunsen Flame 

The data for a mixture of 15 percent propane and 85 percent oxygen 
are according to the foregoing test data and the literature, respective- 
ly, as fellows : 

molal weight: M = 34.1 

density: P a = 1.52 x 10 g cm 2 3 at N.T.P. (=760 mm 3g, and 0° C) 
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Initial concentration: nfc a = 6.7 X 10 ° mol/em 5 fuel E'.T.P. 

flame velocity: wpp = (300 ± 30)cxa/sec) at N ; .T.P. 

initial temperature: d a = 0° C 

terminal temperature: % = 3000° C (hypothetical value) 

ignition temperature : = 500° C and 200° C, respectively 

mean temperature in preparatory zone: -b = 250° C and lOO 1 "' C, 
respectively 

Since it was not possible to make any predictions about the actual 
amount, beforehand, two ignition temperatures were tentatively assumed. 
The specific heats cp were additively computed from the mol heats fir, 

of the single gases. The assumed values are contained in table 4, in 
which the limit of error of the propane values was intentionally put 
very high: 


TABLE 4 


$Q C i 

j 

q^q (reference 9) 

Cn propane * 

l 

1 Cp mixture 

1 j 

1 i 

i i_. 

cal/mol grad 

cal/mcl grad 

cal/mol grad 

i i 

I 100 ! 

7.15 

16.4 ± 3.0 

8.54 * 0.4p 

! 250 ' 

7.50 

24.8 ± 4.0 

10.1 - 0.60 


£ n mixture 

if 

Hence for c_ the numerical values 

34.1 

at d = 100° C: Cp = 0.25 x 0.01 cal g"' 1 grad"* 1 

« = 2p0° C: Cp = 0.30^ O'. 02 cal g — 1 grad -1 

The thermal conductivities of the mixture were computed additively 
from the values for the pure components by means of the following 

•‘•Chemical Jng. vol. Ill, 1 p. 102 gives for 327° C for CH,, C a H 6 
and n— C 5 H 10 the mcl/heats 0^ = 12.6.1— > 20 . 98 -> 38 . 0, hence the value 
29 x 5 for propane at the same temperature is entirely plausible. For 
0 C the propane value was assumed at Cp = lit 3 . The Cp values . 

for temperatures between 0° C and 327° C are linearly interpolated. 


1 6 
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numerical values (X in kcal m * 1 h 1 grd 1 ) , the limit of error being 
again made very wide. 


TABLE 5 


1 

i 

\ o 2 

1 

\ propane 1 

X mixture 

* — 1 

100 

r~ 

0.027 

0.027 ± 0.003 

0.027 ± 0.000 

200 

.031 

.040 + .015 

.032 ± .002 

250 

.033 

.047 ± .020 

.035 ± .003 

500 

.044 + .006 2 * 

0 

0 

# 

0 

1 

.049 ± .010 


Hence for X in c.g.s. units: 


d z = 500° C X z =(136 ± 28) X 10 6 cal cm 1 see" 1 grad 1 

d = 250° X = (97 ± 9) X 10“ 6 

n - 6 

d z = 200 X z = (89 ± 6) X 10 

d = 100° X = (75 ± 2 ) X 10 -6 

•-Frcm equation (10), ( 7 ), (12), and (l4) then follow for the two 
assumed ignition temperatures d z = 500 and 200° C the data of table 6. 


IFor propane, X was linearly interpolated between 0.013 at O'-' C 
and 0.080 at 500° C. The first figure was interpolated from the series 
CH 4 , C 2 H 6 , i— C5H12 with X = 0.026 ->0.0155-^ 0*0105 (Chem. Ing. I, 1, 
1933 , p. 323) and is in agreement with the Landolt-Bbrnstein data. 

The value 0.080 at 500° C was estimated from the kinetic gas formula 
X = ic* 2* Pci, (c*, mean molecular velocity, 2*, mean free wave length, 

Ci, effective specific heat (cf. reference 8, vol III, 2, 1$>25, p. 92) on 
the assumption that at transition from 0° C to 500° C the expression 

1 c*2* P is doubled (exactly as for oxygen, where the effective specific 

heat Cx remains practically constant) and that Cj_ is tripled. 


2 For oxygen, X at 500° C is not known, for nitrogen, it is 0.040. 

The oxygen values are a little higher on the whole; so the value 

0.044 ± 0.006 is probably correct. 
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TABLE 6 


' °c 

&v 

cm 

j ( ^ 

V S ) z 

grad c.m -1 

u 

mol fuel/cm 3 
sec 

Sr 

cm 

1 

i j 

t 1 

500 

t j 

(3.5^0.6)xl0 -3 

: aac. ~ .vor? 7V.. irrizyTr-Ji 1 . 

(5-0^1. 3)xl0 5 

- x(0.i+0 t 0.10) 
f 

fx(5.0±1.5)xl0~ 3 

1 

!200 i 

, ! 

(3-3±0.3)x10 -3 

(2.6i0.3)xl0 5 j 

X(0.l8i0.02) 

n. 

fx(l0.9 £ l.3)xl0 


The cited limits of error correspond to the uncertainty of the em— 
ployed Cp, X and \ z values. 


An objection may be raised against the thus computed values; for 


fir it is questionable whether these Cp and X values, obtained experimen- 
tally on slow processes, are still applicable to the speedier processes 
in the flame. On the one hand, the preheating takes place within about 

3X10' 


— 3 


— 5 


300 


= 10 sec, whereby a single molecule receives only about 10" 


impulses. On the other hand, it is a known fact that a molecule, in 
order to change its vibrational quanta to translation energy must receive 
frem 10 4 to 10 5 impulses, and that this figure is never less except in 
special cases; and that is, when at the impulse, so to say, the complex 
of a chemical compound able to exist is preparatorily formed (reference 
10). It is therefore entirely conceivable that smaller Cp and X ' 
values with only partial or not at all excited vibrational heats are 
valid for the flames. The latter case can be quantitatively computed, 
when it is remembered that, by kinetic gas theory, molal heat Cp and 

thermal conductivity X are represented in the following manner: for 
complete adjustment of the vibrational quanta: 

c p = C t + C r + c 8 + R ( 3.6) 

X = kCj. = k I" 2.5 C t + C r + C s | (17) 

for unadjusted vibrational quanta: 

0} = C t + c r + r (18) 

r~ “1 

X* = kC x = k ; 2.5 C t + C r | (19) 

where K is a constant solely dependent on the mean molecular speed c*, 
the free path length l*, and the density P (reference 8), and R 
is the gas constant Ct the translational portion, C r the rotational 
portion and C s the vibrational portion ®f the molal heat. Quantity C x 
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Identifies the effective molal heat to he inserted in the kinetic gas 
term for thermal conductivity according to Chapman and Enskog. There— 
f ore, in general : 


9 .i - 9 p 

= C ; {~ Cp = 1.5 C t - R = 2.5 cal/mol grad ( 20 ) 

The numerical values obtained by this method are reproduced in table 7 ; 
the values for the gas mixtures were computed additively from the 
values of the individual components by the mixing rule. 


TABLE 7 


c 

Propane 

Oxygen j Mixture: 19 percent prepane] 

85 percent oxygen j 

100 

200 

250 

500 

j 

£1 01 c x C'x 

p — x — x 

16 -4 8.0 18.9 10.5 

22.0 8.0 24.5 10.5 

24.8 8 0 27.3 10.5 

35.0 8.0 37.5 10.5 ! 

: 

^ a s'! ! % si a -cl- 

i | 

| l 

7.15 7.0 9.65 9.5 ; 8.54 7.15 11.04 9.65 

7.38 7.0 9.88 9.5 •' 9. 57 7.15 12,07 9.65 i 

7.5 7-0 10.0 9.5 ! 10.1 7.15 12.59 9.65 

: 8-0 7.0 10.5 9-5 ! 12.1 7.15 14.55 9.65 


Koting that according to equations (iQ)and (l4) 



(21) 


and 



( 22 ) 


and identifying the quantities computed with the unadjusted vibration 
quanta by a dash as before gives for the calculated mixture values of 
table 7 


for ■■«!= 500 0 0: — 8 ^ - C 2a 50 = 1.03 

5T c l 2 so 




and 


^ C Lr.rp- 


= 0.94 




= 200° C: 


o T v 


&v 


5 G 0 P25 G 
^Sso 0 ! 5 00 

C T 1 Ct, 

iioo 21 o o_ -^^04 


(24) 


^100 C lioo 


at * z 


(25) 
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and 


fn - ElSOQ £b2J2Q. = .96 (26) 

R p t n 

E "Pi 00 I200 


Thus the calculated thicknesses of the preparatory and the reaction 
zone vary very little, if in place of the complete vibration adjustment 
only a limited one or none at all is assumed. 


For bR, of course, yet another possibility is involved, namely, 
that in equation (l4) the vibrational heat is completely adjusted for 
Cp tut not for \ z • For according to the physical derivation of equa- 
tion (l4) from equations (7) and (13)> Cp refers to the whole pre- 
paratory zone, and the gas particles remain longer in it than in the im- 
mediate vicinity of the Diace of ign ition where is to be formed. 

For the reaction zone thicknesses 6^ computed with this somewhat ex- 
treme assumption the result is 


at 

at 


$z 


= 500° 


5 R 500 C p250 

C : — * 0* bo 

Sr -Sp 250 — 1 500 


200 ° 



IT 


-13^00 = 0 .80 

•2d .Qi 
iioo i 200 


In this instance the would therefore he smaller than 

viously calculated 03, hut by no means greater. 


(27) 


(28) 


the pre-- 


These considerations show that the zones of preparation and re- 
action defined hy equations (10) and (lk) are still comparatively 
accurately computable, although the data for the specific heats and 
thermal conductivities are not too accurately known, and although it 
is not even decided whether the degrees of freedom of vibration of the 
molecules during the brief intervals in the flame front are excited com- 
pletely, little, or not at all. 


A somewhat closer observation is necessary for the rate of chemi- 
cal reaction U on the one hand, and the thicknesses 5^ and 

on the other. The rate of chemical reaction 

U** 0.2 mol cm sec^ 1 


is quite high; because, when it is remembered that in the reaction zona 
on an average, fuel of the order of magnitude of 2xlC r 6 molecules per 
cubic centimeter is present (at start n^ a = 6.7xl0~ 8 
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2 X 10 ° — S 

a mean combustion period of = 10 sec is obtained. 

0.2 

This period is extremely short compared with tne reaction periods for 
which investigation methods worked out so far, exist. For, in kinetic 
reaction studies by the static method (closed off reaction vessel in 
which the progress of the reaction is followed by observation of the 
pressure variation or some other physical property) time intervals down 
to only about 1 eecond can be followed up; by the conventional flow 
method it succeeds only to about 1(T S or 10~ 3 second; insofar as reason- 
ably defined pressure and temperature conditions in the reaction zone- 
are demanded. Nevertheless, a reaction period of 10 -5 second is still 
not the shortest one conceivable for a combustion process, because a fuel 
particle still receives about 10 5 impulses during this period. Hence, if 

10 2 impulees on a fuel particle were sufficient for its complete decompo- 
sition, which is entirely conceivable with correspondingly energy— rich 
and adequate impulse partners, the rate of combustion would have to be 

10 3 times higher. Such enormously fast combustion processes have, in 
fact, been actually observed, namely, at detonations, where the speed of 
travel of the flame front is counted by kilometers per second; while at 
the normal speeds of ignition in the Bunsen flame a few millimeters per 
second are sufficient. 

Hydrodynamically a comparison of the thickness of the preparatory 
zone &v and the reaction zone bp with the size of the turbulence 
bodies occurring in the turbulent flow (gas bodies with occasionally 
-uniform velocity vector) is interesting. 

Interesting from the point of view of chemical reaction is the com- 
parison of the theoretical guantities 6 V and bp with the experimentally 
directly observed thickness of the luminous zone 5t,. It was established 
at (102) x 10~ 3 centimeter in the previously reported tests, which is, 
within the limits of error, independent of the propane concentration. 

Nor was there any indication that these luminous zones were artificially 
enlarged for the eye or the photographic plate by any flickering motion 
of the flame. In consequence the following possibilities exist accord- 
ing to table 6: 

b T ~ d, r + b^ if* „ = 500° C 
L~VK * I , . 

> assumption A 

and f =1 case I J 

ca se that?> z = 500 ° C 

and f = 2 case II 

^assumption B 

5l & 8^ in case $ z = 200 ° C ' 

= 1 case III; 


and f 
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Which one of these three cases actually obtains has not been definitely 
decided e.s yet, especially as it is not known in what manner the emitted 
spectra are distributed over the individual sections of the luminous 
zone-. 1 2 * But tentatively two assumptions A and B may be made: 

A. Luminosity extends over preparation and reaction zone; in which 
case 500° C should be assumed as ignition temperature, according to 
table 6, that is, a value which, in contrast to previous arguments, ap- 
proximately corresponds to the ignition temperatures obtained by Eernst f s 
method for adiabatic compression or with the Dixon burner. But it would 
be difficult to explain why the luminosity in the preparatory zone, where 
there is practically no chemical conversion, should appear with about the 
same intensity as in the reaction zone. 


B. Luminosity occurs only in the reaction zone: in case II with 
f = 2, v z ~ 500° C would then have to be assumed as ignition temperature 

(cf. table 6), that is, the same high ignition temperature as with assump- 
tion A, or else d z £ 200° C in case III with f = 1, that is, a sub- 

stantially lower ignition temperature • Case II is, however, less probable, 
while it would yield a reaction thickness 5^ which agrees with the 
thickness cf the experimental luminous zone (table 6), this thick- 
ness would not be uniformly luminous, but practically only on the 

first half. For, when it is assumed that luminosity occurs only in the 
reaction zone, its intensity is presumably parallel with the chemical 

conversion and hence with the momentary temperature gradient ail ± n 

ox 

the reaction zone o But its course differs considerably in both 
cases II and III; that is, in case II the high values of the temperature 


gradient occur only on the first half of the reaction zone, while In case 
III the ' temperature gradient is approximately constant over the thickness 
This is readily apparent when the temperature curve for both cases 

is plotted (figs. 7(a) and 7<(b))since absolute value and temperature grad- 
ient for starting point and end point of Sp are known. (See table 6. 

Note that must equal 0 for x = ) 


1 If, for instance, at the beginning of the luminous zone spectra 
other than those of OH, OH, and were found in appreciable quantities, 
an encroachment of the luminous zone upon the preparation zone would be 
quite probable. 


2 0f course, it could be assumed for the present that a particle is 
luminous at a later place than where it is energetically excited, and 
that through it the luminous zone would be wider than the area of the 
maximum temperature gradient. But this type of enlargement is surely 
entirely negligible; because the excited particle has but a lifetime cf 
about 10~ B second, after which it radiates its energy, and in this 
interval it would have traversed exactly l/lOOO of the total reaction 
zone thickness 5p. 
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Comparison of the experimental with the theoretical 5 V and 

8p values leads to the very probable final result: 

The experimentally defined luminous zone 8p either extends to 

the preparatory zone & v and to the reaction zone 8p. In this case 
the ignition temperature in the flame is almost exactly as high 

as by Kerns t } s method of adiabatic compression or with the Dixon burner. 
Or the luminous zone Sp covers only the reaction zone 5p, in which 
instance the ignition temperature i3 z in the flame is lower than for 
adiabatic compression or with the Dixon burner 

Of these two cases the latter is, for the time being, held to be 
more likely. Further elucidation is, first of all, to be expected from 
the spectroscopic side, where type and intensity distribution of the 
spectra in the luminous zone would have to be explored. 


3. Enlargement of a Laminar Free Jet and 
Form of Bunsen Cone Possible in It 


In order to obtain the form of the flame surface of a Bunsen flame 
quantitatively, for laminar flow, the simple free jet without combustion 
is considered and it is proved that, under the previously employed ex- 
perimental conditions (Re > 600), its enlargement at a distance of sev- 
eral tube diameters, corresponding to the height of the measured flame 
cone is very insignificant. 

According to an article on laminar jet expansion by H. Schlichting 
(reference ll) the half— value width 2b 1 of the velocity profile (cf. 
fig. 8, which was taken from that report) depends only on the kinematic 
total momentum K, on the kinematic viscosity v , and on the distance 
x conformably to the equation 


with 



2 . 57 / 



V 


oo 

p 


K = / 2 n u 2 r dr 

J o 


( 29 ) 


( 30 ) 


x In the radial distance b the velocity is occasionally half as 
great as on the jet axis, equal flow cross section being assumed. 
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Of course, equation (29) applies, for the first, only to a jet that 
issues from an infinitely small orifice, hut even so it can be applied 
to the case of a jet from a finite pipe cross section, when the left- 
hand side of equation (29) is regarded as a measure for the jet expan- 
sion. The kinematic momentum which must be the same on all subsequent 
jet cross sections according to equation (30) can then be summarily cal- 
culated for the stationary pipe flow with parabolic velocity profile: 


n 


K = / 2it (2u) ? [ 1 - (g) 2 ] 2 r dr 


= 5, da 


(31) 


By insertion in equation (29) 


-iL. = 2.57 u. = 10 - 28 


(32) 


3 n II s d 2 


Re 


The relative expansion of the laminar free jet issuing at the Reynolds 
number Re from a tube is therefore inversely proportional to it, and 
at the lowest value Re = 612 of the present tests becomes 


2b 10.28 1 

~x _ ~6l2” ~ 60 


(33) 


This is very little, so that there can be no perceptible variation of 
the velocity profile at a distance of a few tube diameters from the 
orifice. The combustion zone will, of course, modify the flow con- 
ditions in the hot part beyond the cone of the flame very substantially. 
But on its inside, that is, between this and the pipe orifice, the 
original velocity profile is scarcely altered; for this would involve 
the variation of a flow pattern before a center of disturbance (that is, 
the surface of the flame), and such changes are usually small, unless 
the center of disturbance necessitates a complete deflection of the flow, 
which, however, is not the case here. 

Thus the assumption that all gas particles enter the surface of the 
flame at the same rate of flow at which they have left the orifice of 
the burner tube appears to be justified. 

Now the shape of the flame surface con be calculated when it is 
assumed that the flame speed w.^ is the same at all points and first 
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1)00013103 abruptly zero at the base. (See reference 4 , p. 46 ; also refer- 
ence 12 .) 


For every single surface element of the flame cone with the horizon- 
tal and vertical coordinates r and h there results in analogy to 
equation (l) : 


Wfi 2n r 


J ■ *£? 


dr = u 2 n r dr 


( 34 ) 


cr, considering the parabolical velocity profile in laminar flow (with 
average flow velocity u) 


u = du 


i J v -> i 

VRy 


( 35 )- 


d(h/R) 


d(r/R) 


Jt j 

\ W 1 \R/ ; 


( 36 ) 


The last relation was graphically integrated for various assumed values 


u 


or — and the cone forms illustrated in figure 9 ascertained. They 

Wfi 

are not in complete agreement with the experimental flame cones, as is 
especially noticeable from figure 10 , where experiment no. 47 — for which 


u 638 , Q 

- - 4.98 and 


wf 1 128 


638 

100.8 


= 6.33 


is exactly valid (cf . table 3) is compared with the theoretical curve 
for — H_ = 6 . Two substantial differences were noted: 

Wfi 


(a) The base of the experimental flame cone not only covers, as in 
the theoretical case, the internal burner tube cross section, but reaches 
beyond it . 

(b) The upper end of the experimental flame cone is, in contrast to 
the theoretical, rounded off and lies considerably lower. 


"’The gas is assumed to be at rest at the pipe wall; whereas 
Michelson, who made the first calculations of this kind, erroneously 
took the possibility of a perceptible sliding into consideration. 
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Finding (a) is probably explainable as follows: immediately abo\e 
the rim of the burner orifice where the ga3, owing to adequate heat re- 
moval by the burner tube, cannot ignite, a considerable radial and 
outwardly directed flow must exist, since there is always a higher 
pressure p a inside of the flame cone than outside (pe) • This radially 
discharging gas cannot ignite until it is a certain distance away from 
the burner tube, and thus the base of the flame cone is drawn cut over 
the- burner orifice. The pressure difference is, according to the mo- 
mentum theorem (reference 4, p. 44): 

Pa ~ Pe = p a w f 2 (p~ ~ ( - 31 ■ 

e 


hence for experiment 47, illustrated in figure 10, 
with 


md 


P a = 1.50x10 3 gram centimeter 3 , wp = 114.4 centimeter sec 


— a 


p a + 273 3273 


= 12 : 


P e b e + 273 273 


Pa - P e i 1 • 50 xlO x 3 1 • 3 lx 10 4 »,11 


dyn 


square centimeters 


= 216 = 0,22 centimeter H 2 0 (S3) 1 

square centimeters 

This pressure difference is quite considerable. If it is regarded as 
dynamic pressure at density P a , a velocity w . can be correlated, which 
should represent .an approximate measure for the radial, outwardly 
directed velocity immediately a.bove tho burner tube rim. Hence from 
equation (37) 


iFor a stoichiometric mixture of l6.7 percent propane, the rest 
oxygen, it would give 

p 

P a = 1.44x10-3 g cm “9 w f = 300 cm s - -, -—£12, and hence 


P a - P e = 1425 = 1.45 cm H 2 0' 

cm 2 
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The last relation ie as jet not to be found in the . literature . Still, it 
is important, because it makes it possible to decide forthwith whether 
or not air from the outside can diffuse in the inside of the cone through 
the free slit between burner tube rim and flame cone base. In general, 
the negative seems to be the case. For on assuming even the very small 
value x = 0.01 centimeter for the radial thickness of the diffusion 
layer and 0 = 0.2 sruare centimeter per second for the diffusion co- 
efficient of air, the inward rat© of diffusion of the air particles 
would be 

= 9 ^ 3 — - go centimeters ner second (40) 

x 0.01 


against the radially outward flow velocity of w = 4.7 :<ll4 = 536 centi- 
meters per. second. 


Naturally, it is also conceivable that the gas portions flowing 
radially outward above the burner tube rim reach combustion only afcer 
appreciable mixing with the surrounding air. They would burn slower or 
faster (slower for initially poor mixtures, faster for initially rich 
mixtures). The recorded flame velocities could he so much more faulty 
as the gas portion, converted in the lowest border zone of the experi- 
mentally surveyed flame cone, is greater. This effect should be watched 
particularly on small burner tube diameters, although it did not seem 
to have any part in the previously reported measurements. The laminar 
flame soeeds obtained with burner tubes l(d = 0.1385 centimeter) and II 
(d = 0.218 centimeter) are satisfactorily grouped within the limit of 
error of about ±10 percent about the dashed curve of figure 5. 

Finding (b), the rounded off and lower position of the upper end of 
the experimental flame cone (cf. fig. 10) is perhaps due to the fact 
that the flame speed wp at the cone tip is higher than on the rest of 
the surface elements of the cone envelope. The gas particles which 
reach combustion in the cone tip receive heat and active chaip carriers 
not only from this tip itself but also from the immediately adjacent 
parts of the enveloping surface, so that ignition and hence combustion 
itself are facilitated. This interpretation is particularly supported 
by the fact that the curvature radii at the upper end of the measured 
flame cones are of the -same order of magnitude as the thicknesses of the 
luminous zones, namely, 0.01 centimeter. 
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4. Prandt 1 1 s Exchange Quantity and Mixing length 
in the Turbulent Free Jet 
(Nuclear and Marginal Turbulence) 

a. General consideration s*— Ax cording to Prandt 1 the developed tur- 
bulence of a stationary flow can be identified in every point of the 
field of flew by two locally dependent quantities, that is, the mixing 
length Z and the turbulent exchange quantity . The former is de- 
fined as length and represents the approximate mean diameter of the tur- 
bulence bodies (that is, areas of closely spaced particles, which oc- 
casionally manifest a uniform vector of the flow velocity); the latter 
is an effective diffusion coefficient of the dimension square centimeters 
per second (reference 13). The introduction of these two quantities is 
based upon the analogy of all transport processes in a gas at rest and 
in a turbulently flowing medium, so far as they are superposed at least in 
the latter case, on the principal flow. These transport processes (that 
is, diffusion = transport of a certain type of particle, viscosity = 
transport of momentum, thermal conduction = transport of energy) are re- 
duced in both cases to irregular motions, in the case of still gas to 
the gas molecule, in the case of turbulent medium to that of the turbu- 
lence body. In still gas the individual gas molecules move irregularly 
at a mean velocity c* and traverse the mean free path length 1* be- 
tween two collisions. In turbulent flow the individual_ turbulence bodies 
possess, aside from the stationary mean flow velocity w, an irregular, 
rapidly changing additional velocity w ! , and they maintain these speeds 
on the average occasionally over a path length Z, the mixing length. 

Just as the diffusion coefficient in the kinetic gas theory, is de- 
fined as 


D ~ i Z* c* 

3 

so Prandtl defines his exchange quantity at 


OD 


e = l w* 


(42) 


| w* | is the time average value (denoted by upper dash) of the absolute 
amount of w f i. The kinetic gas factor l/3 in equation (4l) is included 


1 Averaging over j w* ! is necessary, because the time average value 
of w ! disappears by definition. 
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in equation (42) in the quantities 2 and , which is physically cor- 
rect , since, the mixing length '2 in the turbulent flow must, for 
reasons of continuity, be of the same order of magnitude as the diameter 
of the turbulence bodies, and this diameter cannot be arbitrarily de- 
fined, because every one of these gas bodies is enveloped by a more or 
less thick marginal zone in which the different flow altitudes of the 
adjacent bodies are continuously adopted. 

Numerical values for the exchange quantity c can be secured in 
several ways. For very rough e stim ates it is sufficient, for instance, 
to determine quantity 2 and !w *1 separately in equation (42) (for 
instance, by hot-wire measurements and determination of the correlation 
coefficient for the concurrent velocity fluctuations at two adjacent 
points) 1 and then simply multiply these values with each other. But 
this affords only the order of magnitude of e , Bince the turbulence 
body diameters 2 must be indicated with wider margin. Exact values 
for c are obtainable by surveying a specific transport process in the 
turbulent flow. Only such e values as those obtained by the latter 
method are U3ed in the following study. 

Consider the specif ic__ca8e of a one— dimensional principal flow 
(w = ui + vj + w]c with v = w = 0), superposed by the turbulent veloc- 
ity fluctuations w* = u : i +v*j[ + w'k, and assume that the unit mass ?of 

the flowing medium at any one point has the additive property s — 
expressible by a digit — of varying from place to place. Then, if 

— 0 conformably to figure 11, a "property flow s" transverse to the 

by ~ y 

principal flow direction i can be computed, that is, the amount of 

property s which is transported in unit time per unit surface in the 

j— direct: on 

a = ! v'» ! | Ps - -HP- s ) - ] _ |V“ ps + — ! (43) 

2 by j 2 by _ 

= _ i - iTH . ffi g 2. 

by 

or, with equation (42) taken into considera.tion: 


1 0n a wheat field waving in the wind the turbulence bodies of the 
air current are directly delineated, so that for this case the magnitude 
of the bodies and — by allowance for the mean wind velocity — the 
mean velocity fluctuation) w'j can be given. 
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£y 


(44) 


If, for instarj.ce, s signifies the mean momentum per unit mass in prin- 
cipal flov direction (hence s = u) , then s v 'becomes the ahearit-g 
e tr e 8 -- tr an 8mi 1 1 e d transverse to the principal flow direction, and 

with it for incompressible fluid 

dy 


T = pe 


(49)- 


which means, however, that the exchange quantity e can immediately 
specif led, when the local shearing stress distribution for the particular 
flow and the velocity distribution u have beer, measured. If s in 
in equation (44) indicates the enthalpy s=Cp # 'of a practically incom- 
pressible fluid referred to any one zer^ point. 


•§y ~ ~ c p p € 


§y 


/ \.r \ 

( u oj 


represents a heat flow, and hence yields an effective thermal conductivi- 
ty coefficient in the turbulent flow: 


^turb 


c P £ 
P 


(4?) 


Lastly, if the number of particles of type k, which a" one 

is transported by turbulent mixing motion in a concentration vessel — 


in y— direction per unit time and surface is of interest, then 


s^ r = and Ps = are .written in equation (44), so. that 

■ dnfc ■ 

N k = - e — 

°y 


(48) 


This equation makes it particularly plain that - e is a quantity of the 
type of a diffusion coefficient. 


But in the following not only the exchange quantity e , but also 
the mixing length l is of interest, although it gives the turbulence 
body diameter only in order of magnitude; for l undoubtedly represents 
a characteristic length for the turbulent flow, . just as the preparatory 
zone thickness '5 V and the reaction zone thickness constituted 

characteristic lengths of the stationary flame front in laminar flow. 
Quantity l can be determined from the correlation of the velocity 
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fluctuations on two adjacent points, as previously stressed. However, 
Prandtl and his school used in part a different method. Prandtl had 
found from continuity considerations that in the one-dimensional prin- 
cipal flow 

v* -zifa| (49) 

oy 


must approximately apply; whence for the exchange quantity 

i du 

5F 


5 - 12 l3y 


(50) 


or by insertion in equation (45) : 

T „ p Z 2 | ^ | 

- Oy oy 

In the last equation, t, p and — are experimentally accessible, sc 

Sy 


(51) 


that Z can be computed immediately. 

Equations (45) and (5l) formed the initial equations by which the 
quantities e and Z, in the subsequently discussed reports, were 
determined. 


(3. Nuclear and marginal turbulence in a free .let .— The interpreta- 
tion of the flame phenomena observed in turbulent flow stipulates e 
and Z as functions of the location for the turbulent free jet issuing 
stationarily from a tube. Although this case has up to new not been 
explored experimentally in its entirety, two other cases are known: 


a. The turbulence in a circular pipe in stationary flow was ex- 
perimentally investigated by J. Nikuradse and the e and Z values 
determined (reference 14) . 

b. The turbulence in the mixing zone between a turbulent free jet 
and the bordering still air was theoretically explored by Tollmien 
(reference 2); he found the sole constant still undetermined in the 
theory by comparison with Gottingen wind-tunnel measurements (reference 
15), so that here also, the two factors Z and e , are known. Thus 
the turbulence conditions in the free jet can be immediately surveyed 
when assuming that the jet can be decomposed in a nuclear and a marginal 
flow, and that in the former the unchanged turbulence condition of the 
simple pipe flow is still maintained (according to a). This decompo- 
sition cannot, of course, apply to any desired great distance x from 
the tube orifice, because the marginal zone expands continuously, while 
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the nuclear flow shrinks conically. But it will he seen that an essen- 
tial part of the turbulent combustion phenomena (and this is precisely 
the most interesting part) still falls in the zone of the nuclear flow. 


Figures 12 and 13 show the nuclear turbulence in a free jet in dia- 
grammatic representation. In figure 12 the dimensionless exchange quan- 
tity — (u is flow velocity averaged over the tube cross section and 
uR 

with respect to time; R, tube radius) is plotted against the relative 

distance from the tube axis -; 1 in figure 13 the mixing length |j- 

E *> 

r 

made nondimens ional with the tube radius is alee plotted against — . It 

R 

is seen from figure 12 that € is not maximum at the tube center line 


but approximately at distance r = ^ R. This fact can be regarded as 


an experimental confirmation of the approximate formula (49); for the 
closer the approach to the tube axis, proceeding from the wall, the 

<3u 


greater the reduction in 
e muet again decrease. 


'dy 


and since l cannot grow indefinitely. 


According to figure 13 the mixing length is less than 8 percent of 
the tube diameter, hence less than 0.l6 millimeter for a 2-millimeter 
burner tube. Accordingly the turbulence bodies would, in consideration 
of the nonarbitrary definition of their diameters, definitely come 
within the order of magnitude which had been found for the thickness of 
the laminar flame zones. 

The nondimens ional quantities y- and /y- employed later, 2 are re- 
produced in table 8 for nuclear flow at different Reynolds numbers and 

different — , 

R 


1 Nikuradse , s graphical representation could not be used direct, be- 
cause he used the so-called "shearing stress velocity v*" to make e 
dimensionless, which, however, is not so suitable for the present purpose 
as the mean flow velocity u. 

2 According to the kinetic gas theory the kinematic viscosity y 
is approximately equal to the diffusion coefficient D; for a pure gas 
the agreement is complete. 
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The marginal turbulence at the boundary of a free jet van explored 
by Tollmien on the streamline pattern illustrated in figure 14. This 
does not completely correspond to the boundary conditions of the present 
tests, which involved the exit of a free jet from a tube of finite wall 
thickness and tapered at the orifice, rather than the discharge of a free 
jet from a straight wall (visualized as being placed perpendicular to the 
plane of drawing through the y-axis) as in figure 14. Nevertheless, an 
analysis of Tollmien *s case is interesting. The two-dimensional plane 
problem which Tollmien first investigated and which presents tho greatest 
interest, involved the solution of the differential equations 


and 


Su bv - 

bv. ou _ 1 ^ T xy 

u — + v — - — ; — 

ox by P by 


with the shearing stress 


T xy = Pk 2 x 2 ^ 

7 Sy Sy 


(52) 

(53) 


(54) 


Herein ' k is a temporarily unknown constant in the formula for the mix- 
ing length 


l = k x (55) 

The differential equations were solved by means of a stream function for 
the velocities u and v, without going into any further details. Only 
the most important result is repeated here: namely, the entire marginal 
zone lies in a sector, the tip of which is located at the tube end (in- 
side wall of tube) and this is bounded by the two straight lines 


3 



Jx = 0.981 ✓ 2k 2 x = ti 

(56) 


3 

y 2 = — 2.04 1 / 2k 2 x = — b a 

(57) 

For 

the 

y > yi the flow velocity is u = 5, for y < y 2 , u = 0. 
total width of the mixing zone distant x from the mouth 

Hence for 
of the tube 


3 

b = b x +b a = 3-02 / 2k 2 x 

(56) 


The best agreement with the experiments of the Gottingen wind-tunnel meas- 
urements was obtained with 
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^2k 2 =0.845 and k = 0.0174 (59) 

which gives a mixing length 

l = 0.0174 x = 0.0682 t (60) 

The mixing length in the marginal zone therefore amounts to about 7 per- 
cent of its momentary thickness, while Z in the nuclear zone of the 
free jet was established at less than 8 percent of the tube diameter. 

In addition 


bi = 0.83x; b a = 0.172x; b = 0.255x 


(6l) 


From the first equation (6l) it follows that for x = 12 R = 6 d the imar 
part b x of the marginal zone is equal to R, that is, the marginal zone 
here reaches the original tube axis, so that the area of the nuclear flow 
with the original turbulence conditions in the tube disappears. 


Inside of the mixing zone b the flow velocity u drops in first 
approximation from its original value u in nuclear flow to the value 
0 in still air, thus realizing the trapezoidal velocity profiles for the 
free jet, as experimentally observed. With this the turbulent exchange 
quantity is 


P ou ~ Z 2 u 0.0174 2 UX'. 3 = 

= 1 5y b~ = 07255 ~ ~ 1 - lc ® x 10 u x 


(62) 


or nondimens ionally expressed, when the Reynolds number of the tube flow 
R = is introduced: 


= 1.188 x 10 3 X Re - 
u d 


( 63 ) 


This treatment of the free-jet marginal zone is, of course, not com- 
pletely exact, inasmuch as the calculation was based upon the plane case 
in contrast to the really interesting marginal zone which is curved at 
right angles to the x, y— plane of figure 1^, Eence, the closer the 
inner part b^ of the marginal zone approaches the radius of the tube^ 
the greater the discrepancies between the computed e and l values 

that must be expected. Owing to this very uncertainty the exact \.^~ 

i py 

values were not written in equation (62), although it would have been 
entirely possible on the basis of Tollmien^ results, tut the average 


value 


— \~ — . This also meant a great saving in paper work, 
i Py I b 


On the 
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other hand, the exchange quantities and mixing lengths calculated by 
equations ( 62 ) and ( 60 ) cannot he too far wrong, for in the report cited 
Tollmien also treated the case of the rotationaily symmetrical free, jet 
that emerges from an infinitely small orifice. The mixing zcne then 
assumes the shape of a cone with straight enveloping eurface. The mixing 
length at distance x from the point of discharge is 

Z r = 0.0138 x (64) 


and the momentary cone radius 


r 1 = 0.2l4 x 


(6 


0 ) 


This mixing length differs from that used above (Z = 0.0174 x) by only 
9 percent, while the cone radius r* differs by a mere l6 percent from 
the previous width of marginal zone b = 0.255 x. 1 And the exchange 
quantities computed with the values (64) and ( 65 ) differ even less from 
the others, that is, only 2 percent: for the insertion in equation (62) 
gives 


Z»2u 0.0158 2 ux , rp -3 = 

= 1.168 xlO u x 


0.214 


( 66 ) 


i € 

The values of |T and ^ computed by equations (60) and ( 63 ) are 
givey in table 9 for different Beynolds numbers cf tube flow and for 
different relative distances — from the mouth of the tube. Column 2, 

hi R 

containing — , shows immediately the extent to which the marginal zone 
B 

has penetrated into the area of the original nuclear flow. With the 
data of the nuclear zone on top, and the data on the marginal zone 
below, the turbulence properties in the two-flow regions can be directly 
compared. An even better insight is afforded from figure 15, where the 
nondimens ional exchange quantity 4 (solid curve) and the nondixaeneional 

mixing length ~ (dashed curve) for nuclear and marginal f low are shown 

plotted against different — . 2 It is interesting to note that in the 

R x 

vicinity of the tube orifice (for — < 4) the turbulence in the nuclear 

R 


1 TRe cone radius r r must be compared with the thickness b for 
the reason that at the place where the region of the nuclear flow is ab- 
sorbed by that of the marginal flow, the total width of the free jet is 
exactly 2 r 1 and 2 b, respectively. 

2 

The discontinuous transitions between marginal and nuclear zone 
and between marginal zone and the outer quiescent medium manifested by 
the e and 7 curves in figure 15 are attributable to the simplifying 
assumptions in the present estimation, which naturally describes the 
true conditions only in first approximation- 
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zone is greater than in the marginal zone, loth as regards exchange quan- 
tity e and mixing length Z. But as the distance from the tube orifice 
increaees, the € values of the marginal zone very quickly overtake 
those of the nuclear zone, that is, the more the latter is absorbed by 
the first, the closer its mixing lengths approach. At the point where 
the nuclear zone disappears the mixing lengths in the nuclear and mar- 
ginal zone are practically alike. This fact is noteworthy, for it means 
that the numerical values derived from two entirely independent reports 
by Tollmien and Nikuradse are not merely compatible but even supplement 
each other. This is a necessary condition if the turbulence estimate 
which is based on the estimates in two partial zones is to represent the 
actual conditions to some extent. 


5 . The Several Theoretically Conceivable Turbulence Effects on 
Flame Velocity and Interpretation of the Experimental Test Data. 

To understand the flame phenomena in a medium of turbulent flow it 
is necessary to proceed from the fact that the turbulent flow contains 
individual gas bodies which execute an irregular vibrational motion and 
with it enter the combustion zone. The approximate diameter of the gas 
bodies is given by the mixing length Z, the approximate vibratory mo- 
tion by the mean velocity fluctuation |'w’ |. The transport processes 
amplified in the turbulent flow are identified by Prandtl’s exchange 
quantity 

e = Z i s? I ( 67 ) 

defined as a diffusion coefficient. 

It is easy to visualize two turbulent flows 1 and 2 which, in 
macroscopic transport processes (that is such passing over the length 
of several diameters of turbulence bodies) are identically effective; 
hence 

e 1= e 2 (68) 

but at the same time can differ substantially in the hydrodynamic struc- 

ture of the two flows. For condition ( 68 ) does not as yet indicate \ 
that it be necessary that 

li = I 2 ( 69 ) 



and 


( 70 ) 
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Merely the equality of the product l x |w£ I = 1-2 \m.s ! is demanded. 

This indicates that the turbulence of a flow must be identified by two 
numerical values. Previously these were the Brandt 1 quantities . I and 
e; now they are the quantities l and | jjfi . It is therefore necessary 
to distinguish between a "coarse— body" turbulence with small velocity 
fluctuations and a "fine— body" turbulence with great velocity fluctua- 
tions. The effect of these two types of turbulence on the flame phe- 
nomena may be an entirely different one, although under practical test 
conditions both orders of magnitude of turbulence bodies are likely to 
be always present. 

The effect of coarse— body turbulence with l » St- i on a combus- 
tion zone can be explained as follows: Consider a pure laminar flow 

with a velocity profile as in figure l6(a). In this flow a perfectly 
smooth flame surface can exist at a place when the flow velocity u 
and the laminar flame velocity wpp are exactly the same. Then 
visualize in the flow before entering the zone of the flame, the 
sudden appearance of a few turbulence bodies, so that for a moment the 
velocity profile shown in figure l6(b) is created. The flame sur- 
face at section a is now subjected to bulging in the flow direction 
similar to that on the Bunsen burner, where the laminar flame cone be- 
comes so much mere pointed as the flow velocity (u + u* in figure 3.6(b)) 
at point a exceeds the normal flame velocity wfp. But in section b 
where the gas enters the flame surface with a velocity inferior to wfp, 
it speeds against the gas (as in the flareback of a Bunsen burner). An 
inverted flame cone begins to form which grows continuously and thus 
speeds ever faster toward the gas flow. A stoppage or reversal of this 
flame motion occurs only with the sudden appearance of gas bodies at con- 
siderably greater flow velocity in section b, which must substantially 
exceed both u and u + u' , as assmued on point a. Thus the coarse— body 
turbulence roughens the perfectly even flame surface in laminar flow and 
may leave it ruffled. But at the same time the effective flame surface 
apportioned to a definite flow cross section q in the turbulent flow 
becomes substantially greater than in laminar flow, as is readily apparent 
from a comparison of figures l6(b) and l6(a). In consequence the rough- 
ened or fringed flame surface can be kept stationary by a turbulent flow 
at a point in space only when at that point the mean flow velocity u, 
which, indeed must be equal to the turbulent flame velocity wft. referred 
to mean— flow cross section, is greater than the laminar flame velocity 
wfi. 


1 5p is the experimentally determinable luminous zone thickness in 

the flame and indicates the thickness of the laminar combustion zone in 
order of magnitude. 


NACA TM No. 1112 


37 


Thus the coarse— body turbulence effects a rise of the flame .veloc- 
ity referred to flow cross section, even if the flame velocity in the 
individual microscopic flame surface elements, which for the most are 
diagonal to the principal flow direction, is entirely the same as in 
laminar flow; namely, wf^.- 

This concludes, the qualitative discussion of the effect of coarse— 
body turbulence. As to quantitative data it would entail the calcula- 
tion of the extent to which the originally flat, laminar flame surface 
(cf . fig. 1 6 ) is enlarged by the sudden appearance of velocity fluctua- 
tions of the order u l in the flow and when the turbulence bodies have 
on an average the diameter Z. Such a calculation would, of course, 
be almost impossible in the case of completely ruffled flame surface. 

But as long as the surface is only slightly wavy or roughened, such as 
is represented in figure l 6 (b), the flame surf axe could be approximated 
by a series of laminar Bunsen cones disposed side by side wi t h v ar o ou sly 
directed cone tips, similar to the Maker burner where the entire sur- 
face of the flame consists of a series of side-by— side small Bunsen 
cones but all arranged in the same direction. Just as the enveloping 
surface of the flame cone on the Bunsen burner is in first approximation 
proportional to the mean flow velocity in the burner tube, so also it 
must be assumed that the surfaces of the individual positive and negative 
partial cones on the roughened flame surface are proportional to the ve- 
locity fluctuations ±u f * Averaging over the entire flow cross section 
then would yield the proportionality 

Wft ~ I u«! ( 71 ) 

or for constant turbulence-tody diameters Z, which in general depend 
cnly on the vessel’s boundary (cf. fig. 13 , which applies to turbulent 
tube flow) : 

Wf t ~ e ( 72 ) 

£ 

or finally, if is approximately constant (as also holds true 

VRe 

approximately for tube flow = nuclear flow in free Jet, fig. 12) : 

Wft ~ Re ( 73 ) 

On these assumptions, about the occurrence of which in arbitrarily 
shaped reaction chambers nothing can, of course, be stated as yet, the 
turbulent flame velocity should be proportional to the Reynolds number Re. 

What does the experiment say to this? Among the previously de- 
scribed experiments the coarse-body turbulence undoubtedly plays the 
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earliest part in the larger burner tubes II and III with their maximum 
mixing length computed at l = 0.17 and l = 0.21 millimeter for tube 
diameter d = 2.l8 and 2.718 millimeters according to figure (13), while 
the thickness of the experimental luminous zone in laminar flow 

amounted only to about half, that is, 0.10 millimeter. Considering the 
maximum flame speeds obtained for tubes II and III, they manifest 

no direct proportionality with Re, but still an approximate linear re— 
lationehip (fig. 6) 


w fmax = A X Re + B (74) 

The first term on the right-hand side possibly comprises the effect of 
the coarse-body turbulence. But in addition to that there is a second 
effect in term B, which also raises the flame speed and which carries 
particular weight at the small Reynolds numbers such as occur on the 
narrowest tube. Here, however, the turbulence bodies were already of 
the same order of magnitude as the thickness of the luminous zone &p; 
the maximum mixing length computed for tube diameter d = 1.385 milli- 
meters was ?. = 0.11 millimeter according to figure (13). Hence it is 
very likely that the "fine— body" turbulence, at which the turbulence 
bodies are smaller than the thickness of the laminar 6^, also played 

a part in these tests. 


In analyzing the effect of this type of turbulence in the extreme 
limiting case Z« no further roughening of flame surface may be 


expected, because the turbulence bodies have much too small a diameter 
for it. But a different effect must occur, namely, an amplification 
of every microscopic transport process in the flame surface, especially 
between the reaction zone and the preparatory zone. These transport 
processes, which, for example, effect the transfer of active particles 
or of heat from the reaction zone to the preparation zone, govern the 
flame speed very decisively, as is. readily apparent from equation (15)* 
According to it the flame^speed is proportional to the square root of 
the transport quantity — — ; hence 

c p p a 


w 


r fl 


c-oPa 


(75) 


Up to a factor departing little from unity this quantity ~ 


c p p a 


for gases, equal to the diffusion coefficient D and the kinematic vis- 
cosity v . Observing, further, that in turbulent flow, according to 
equation (47) , this transport quantity should be replaced by the Prandtl 
exchange quantity e, it follows that in the extreme limiting case of 
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fine— body turbulence the ratio of turbulent to laminar flame speed is 

*11 = (76) 

Vf p V 

This formula is suitable for a mathematical check of the present study, 
since it not only presents a proportionality like equation (73) (which 
represents a limiting case of the coarse— body turbulence ) , but also a 
real, numerically assessable equation. In this instance the conditions 
are much better for the fine-body turbulence than for the coarse— body 
turbulence. 


If this concept that for small burner tube diameters the fine-body 
turbulence also played an essential part along with the coarse-body tur- 
bulence, is correct, the maximum flame speeds w^ Eax observed in the tur- 


bulent tests should not depart abnormally, in order of magnitude, from 
the flame speeds obtained by equation (76) and from the cited ^ values 
in the free jet. This calculation is made easier by the fact that only 
the (y values of the free Jet nucleus need to be taken Into account, 
which by assumption were to agree with the values in pipe flow; for the 
greatest flame speeds in the present tests are undoubtedly assumed on 
the inside boundary of the turbulent flame zone, and this inner boundary 
definitely lies in all tests (except tests 52 and h5) within the nuclear 
zone of the free jet, as indicated by the dashed line3 of figure 15. 
According to table 3 the ratio of the corresponding experimental inner 
cone surface lines to burner tube diameter, that is to say, was con- 


sistently much below 6. Hence an inside boundary of the turbulent com- 
bustion zone can be calculated by the differential equation 

\ 2 


d(r/f?) 


1/ V 


5 

] „_u ••_\,bma2L 


- w. 


w fl 


J 


£ 

V 


- 1 


( 77 ) 


This equation is in perfect accord with equation (36) except for 

the difference that the values fcr turbulent flow replace the 

u max 

parabolic velocity distribution in laminar flow, and turbulent flame 
speed w fl of equation (76) substitutes for the laminar flame speed 

w ft* Equation ( 77 ) was graphically integrated for such — 

values and those Reynolds numbers which approximately corresponded 
to the foregoing tests and for which the £ values could be taken 
directly from Nikuradse's report (reference it) . The numerical values 

of Re and employed are contained in table 9 - The flame cones 

obtained by integrating ( 77 ) are reproduced in figure 17. As for the 
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evaluation of the experimental cones the length. S. of the enveloping 
surfaces and their center of gravity radius E* (cf. the figs, in the 
4 th and pth cols, of table 9) were determined and from these data the 
flame speed wpt averaged over the entire cone, was easily computed, 

because it is necessary that 


Wft = 


uE g n 
2it S E ' 


(78) 


or, when .formed, in ratio to the laminar flame speed: 

“ d2 


W. 


ft 


R" u 
X 


wfi 2 R 1 S Vfi 


(79) 


These data are reproduced in the next to the last column of table 9. 

They manifest no systematic motion with — — at constant Eeynolds 

wfi 

numbers, so that averages are permissible (cf. last col., fig. 9) . 


TABLE 9 


1 j 

Ee 4xl0 3 

« j — f 

i ! ! 

6.1X10 3 i 16.7X10 3 1 23,3X10® 

Ae j 63.2 • . 

•78.0 ! 129 .I j 152.7 

u/w fl ; 10 15 

10 15 i 20 40 

*{ i 

20 40 

! 

s/E 1 3. 925 5.925 

! 

3.265 4.910 | 4,310 8.59 

. . _ .... _ 1 ... .. . ... 1 

! 3.61 7.41 j 

E'/B 10.435 0.440 

_ i 

o.46o 0.450 j o.46o 0.450 j o.46o o.46o ! 

i i j 

Vft = __S. JL. ;2.93 2.88 

Wfi 2S E' W-p-, j 

3.33 3.39 i 5-04 5.17 i 6.03 5.87 

i 1 

j ; 

i 

j 

[Average for ! 2.90 

1 "n! 

L L 

1 i 

! 

3.36 j 5.10 

! 

1 

} 

1 

j 

1 

5.95 

1 


1 The departure cf the individual values hf_k from the av^rases is 

w fl 

even less than 2 percent, hence lies definitely within the range of 
errors of the employed method of graphical evaluation. 
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The turbulent flame speed, thus computed therefore represents quan- 
tities computable from the entire surface of a flams cone located in the 
nuclear zone of the turbulent free jet and on which the flame speed of 
every flame surface element is increased in the ratio relative 

to the laminar flame speed. Flame speed Wft should agree then in 
order of magnitude with the experimental values of wf max , if a sub- 
stantial effect on the flame speed were due to the fine— body turbulence. 
And this actually appears to he the case: a comparison of the 


w fmax 

values (dashed curve, fig. 6) with the experimental values ~ — ~ 

w fmin 

shows the latter values to differ only by 30 percent from the former 
values; hence the order of magnitude of both quotients is definitely 
the same. Complete agreement is, of course, not to be expected; for 
the assumption of the extremely fine— body turbulence, where the tur- 
bulence bodies should be substantially inferior to the thickness of 
the laminar flame surface, is not fully complied with in the previously 
reported tests. _ 


Wft 


Figure 6(b) further shows that — ~ (/ Re, as established in 

Wfi w fmax 

first approximation for the experimentally defined — values 

wfmin 


at small Reynolds numbers. Admittedly, the factor of proportianality 
for both quotients is different, which proves again that another effect 
besides ' the fine— body turbulence, namely, the coarse— body turbu- 
lence, was involved in the tests. This is borne out by the fact that 
the form of computed inner boundary cones, illustrated in figure 17, 
does not exactly agree with the experimental, with their more bulging 
cones, as is evident from the flame photographs of figure 4. 


V -f X£3X 

Another unusual fact is that the observed — values are smaller 

wft wfmin 

than the theoretically calculated values. But this means that in 

Wfl 

transport processes over path lengths equal to or less than Prandtl's 
mixing length 2, the exchange quantity e - is only partially effective. 
This finding may now seem perhaps trivial, but up to now the author has 
run across no experiment from which such could have been accurately 
deduced. 


In conclusion, a few words about the outer demarcation of the tur- 
bulent flame zones. They extend in their upper part, especially on the 
largest of the employed tubes (ill) far beyond the inner cross sec- 
tion of the tube. Here also the pressure difference before and after 
combustion may be involved as in the enlargement of the base of a 
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laminar Bunsen cone. However, some other factor Besides this potential 
pressure effect must also have been effective, and this is very likely 
to be found in the marginal turbulence of the free jet. For the outer 
boundary surfaces of the turbulent flame zones lie to a great extent 
in the area of the marginal turbulence in the tests with tubes II and 
III. This is readily apparent when comparing the forms of these sur- 
faces of demarcation with the cone— shape nuclear zone of the free jet 
of figure 15, and, further, when reflecting that the ratio hi for 

these boundary surfaces, on changing from tube I~i> tube II tube III, 
moves continually nearer to the value 6 or even exceeds it (cf. 
table 3)* The numeral 6, however, describes the ratio of enveloping 
surface to tube diameter (cf. fig. 15, also table 8) for the cone- 
shape nuclear zone of the turbulent free jet. In tests 52 and 45, 

where - = 12.5 and 7-3* the ccneo extended especially far into the 
d 

area of marginal turbulence. 

Returning to the original question — that is, of the effect of the 
turbulence of a flow on the flame speed — the answer can be worded as 
follows: the turbulence always increases the flame speed referred to 

the average flow cross section, and in twofold manner: for coarse— 
turbulence bodies exceeding the laminar flame thickness, the flame sur- 
face is effectively enlarged, which is, by roughening or ultimate 
ruffling. For turbulence bodies smaller than the laminar flame thick- 
ness this phenomenon probably recedes in the background, though the 
flame speed increases by reason of the fact that all transport processes 
between reaction zone and preparatory zone in the flame surface are 
amplified, owing to the turbulent mixing motion. The effect of this 
fine-body turbulence can now be determined beforehand, but not that of 
the coarse— body turbulence. In the engine the coarse— body turbulence ' 
is most likely to assume the controlling part. Fine-body turbulence 
can be achieved only in flows issuing from very narrow tubes or slits 
of less than 1 millimeter, and such cross sections are avoided as much 
as possible for high— flow velocity because of the substantial flow re- 
sistances incidental to it. 


Translation by J. Vanier, 
National Advisory Committee 
for Aeronautics. 
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Table 2. 


Test 

No . 

C,H. | 

Flow velocity 

| 0, | C,H, + 0, - G 

N . T.R cm'/s 

% C,H, 

Jp 

mm H,0 

G 

dp 

1 - 

2 , 122 - 10 -*- ** 
g cm -1 b~ 1 

9 

g/cm* 

V 

cm*/s 

1 

0,377 

1,949 

2,326 

16,2 

203,7 

1 , 14 - 10 -* 

1 , 86 * 10- 4 

1 , 43 * 10- 3 

0,130 

2 

0,659 

1,475 

2,134 

30,9 

167,9 

1 , 27 * 10 -2 

1 , 67 * 10- 4 

i ; 52 * io - 3 

0,110 

3 

0,862 

1,154 

2,016 

42,8 

147,2 

1 , 37 * 10~ 2 

1 , 55 * 10- 4 

1 , 59 * 10- 3 

0,098 

4 

1,269 

0,754 

2,023 

62,7 

128,1 

1 , 58 * 10- 2 

1 , 34 * 10 -4 

1 , 70 * 10 -3 

0,079 

5 

1,513 

0,516 

2,029 

74,6 

119,0 

1,71 * 10- 2 

1 , 24 * 10- 4 

1 . 77 - 10- 3 

0,070 

6 

1,897 



1,897 

100,0 

97,7 

1,94 JO ' 2 

1 , 09 * 10- 4 

1 , 92 * 10- 3 

0,057 

7 


1,965 

1,965 

0 

188,1 

1 , 04 * 10" 2 

( 2 , 05 * 10 - 4 ) 

( 1 , 34 * 10 - 3 ) 

( 0 , 153 ) 

8 



1,656 

1,656 

0 

158,5 

1 , 04 * 10 -2 

( 2 , 05 * 10 - 4 ) 

( l , 34 - 10 - 3 ) 

( 0 , 153 ) 

9 

0,174 

1,888 

2,062 

8,4 

191,0 

1 , 08 - 10 -2 

1 , 96 * 10- 4 

1 , 39 * 10- 3 

0,141 

10 

1,412 

0,383 

1,795 

78,7 

103,0 

1 , 74 * 10- 2 

1 , 22 * 10' 4 

1 , 80 * 10- 3 

0,068 

11 



1,658 

1,658 

0 

164,0 

1,01 *io - 2 

( 2 , 05 * 10 - 4 ) 

( 1 , 34 * 10 - 3 ) 

( 0 , 153 ) 

12 

1,014 

1,060 

2,074 

48,9 

145,5 

1 , 42 - 10- 2 

1 , 49 * 10- 4 

1 , 62 * 10- 3 

0,092 

13 

1,913 

0,212 

2,125 

90,0 

114,8 

1 , 84 * 10- 2 

1 , 15 * 10“ 4 

1 , 86 * 10-3 

0,062 

14 

1,932 



1,932 

100 

99,5 

1 , 94 * 10- 2 

1 , 09 * 10 -4 

1 , 92 * 10-3 

0,057 

15 


2,000 

2,000 

0 

191,0 

1 , 0410- 2 

( 2 , 05 * 10 - 4 ) 

( 1 , 34 * 10 - 3 ) 

( 0 , 153 ) 


80 * 


Table 8 


Nuclear flow 



Re =4,0-10* 

Re = 6,1* 10* 

Re = 16,7- 10* 

Re = 23,3-10* 

r 

~R 

l 

£ 

1/.L 

l 

f. 

ill . 

l 

e 

i/z 

l 

£ 

i/- 


R 

V 


R 

V 

\ *• 

R 

y 

r •• 

R 

V 

\~ 

0,0 

0,161 

( 2 , 0 ) 

1,41 

0,157 

( 2 , 0 ) 

1,41 

0,154 

( 2 , 0 ) 

1,41 

0,148 

( 2 , 0 ) 

1,41 

0,1 

0,159 

7,1 

2,66 

0,156 

10,1 

3,17 

0,145 

22,1 

4,70 

0,144 

29,4 

5,42 

0,2 

0,152 

9,9 

3,14 

0,152 

13,8 

3,72 

0,143 

30,9 

5,56 

0,143 

41,2 

6,42 

0,3 

0,149 

11,5 

3,40 

0,145 

16,2 

4,03 

0.140 

37,0 

6,08 

0,137 

48,6 

6,97 

0,4 

0,142 

12,7 

3,56 

0,137 

17,9 

4,23 

0,130 

39,8 

6,31 

0,129 

52,6 

7,25 

0,5 

0,132 

13,1 

3,63 

0,130 

18,7 

4,32 

0,122 

42,2 

6,50 

0,122 

55,5 

7,45 

0,6 

0,118 

12,9 

3,59 

0,116 

18,2 

4,27 

0,112 

41,9 

6,47 

0,109 

54,5 

7,38 

0,7 

0,100 

11,8 

3,43 

0,098 

16,7 

4,08 

0,096 

38,8 

6,23 

0,092 

49,7 

7,05 

0,8 

0,076 

9,6 

3,09 

0,075 

13,7 

3,70 

0,072 

31,2 

5,59 

0,070 

40,4 

6,35 

0,9 

0,046 

6,2 

2,48 

1,00 

0,045 

8,6 

2,94 

0,043 

19,8 

4,45 

0,043 

26,2 

5,12 

1,0 

0,000 

( 1 , 0 ) 

0,000 

( 1 , 0 ) 

1,00 

0,000 

( 1 , 0 ) 

1,00 

0,000 

( 1 , 0 ) 

1,00 


Marginol flow 



For any Re 

Re = 4.0-10* 

Re = 6,1- 10* 

Re = 16,7-10* 

Re = 23,3-10* 

R 

b i 

ba 

b i + b a 

l 

£ 

ill . 

£ 

iii_ 

£ 

V * 

£ 

i /* 


R 

R 

R 

R 

V 

\ V 

V 


V 

r " 

V 

\T 

2 

0,17 

0,35 

0,51 

0,035 

4,8 

2,17 

7,3 

2,69 

19,9 

4,46 

27,7 

5,26 

4 

0,33 

0,69 

1,02 

0,070 

9,5 

3,08 

14,5 

3,81 

39,7 

6,30 

55,4 

7,44 

6 

0,50 

1,04 

1,54 

0,104 

14,3 

3,78 

21,8 

4,66 

59,6 

7,72 

83,1 

9,1 

8 

0,66 

1,38 

2,04 

0,139 

19,0 

4,35 

29,0 

5,39 

79,4 

8,91 

111 

10,5 

10 

0,83 

1,73 

2,56 

0,174 

23,8 

4,87 

36,2 

6,02 

99 

10 

138 

11,7 

12,08 

1,00 

2,08 

3.08 

0,210 

28,7 

5,35 

43,8 

6,61 

120 

11 

167 

12,9 
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Table 3 



Mole. 

weight 

it °F >• 
propone 

Flow velocity 
cm % fg.T.P. 



Goging of flame zone 


Lineor 
ve locity 

Flame 

velocity 

max 

Reynolds 

number 


Flame 


Test 

No. 




% 

Propone 






VRe” 




Propane” 

0 2 

G= Pro- 
pone +0^ 

C 

Method*) 

S-Cone sur- 
face line 

R'= c.p. 
radius of S 

F = cone 
Surface 

cm a 

cm/s 

Wp = iL 
cm/s 

^fmin 

Re=^ 

Noisless 

Whistles 

o 

d 









Tube i: (l = 

0,1385 cm 









11 

45.7 

2,12 

21,54 

23,00 

8,90 

0,140 

K 

0,405 

0.040 

0,117 

1571 

202,0 


1 555 

39,4 

yes 

no 

3.30 

12 


2,50 

18,24 

20,83 

12,43 

0,135 

K 

0.312 

0,040 

0.078 

1383 

200,0 


1421 

37,7 

n 

ii 

2,25 

22 

40,5 

7.10 

50,85 

57,95 

12.25 

0,135 

K 

0,075 

0,055 

0,233 

3860 

248,5 


3900 

02.9 

no 

yes 

4,87 

25 


5,08 

38,20 

43.28 

11,72 

0,130 

K 

0,005 

0.041 

0,171 

2874 

252.7 


2933 

54,1 

•i 

n 

4.80 

49 

40.17 

2.028 

13.05 

1 5.08 

10.75 

0,120 

Ph 

0.245 

0,243 

0.037 

0.032 

0,057 

0,0488 

1042 

275.0 

321.0 

1,17 

1119 

33.5 

yes 

no 

1,77 

1.75 

50 


3.04 

10,33 

13.97 

20.02 

0.117 

Ph 

0,300 

0,350 

0.0425 

(0.031) 

0,090 

(0,068) 

028 

145,2 

(205.0) 

(1.41) 

1098 

33.1 

• i 

• i 

2,00 

2,53 

51 


7.00 

29,55 

37,45 

21.10 

0.123 

Ph 

0,565 

0,538 

0.050 

0,038 

0,177 

0,128 

2488 

211,4 

291,9 

1.38 

28(H) 

52,9 

no 

yes 

4,08 

3,88 

52 


10.34 

27,85 

38,19 

27,07 

0,110 

Ph 

1.73 

0.07 

0,70 

2537 

50.2 

V 

3035 

55,1 

•I 

1 1 

12,5 

53 

40,10 

7.59 

45,00 

52,59 

14,45 

0,132 

Ph 

0,500 

0,438 

0.057 

0.044 

0.200 

0.121 

3500 

202,0 

434,7 

1,00 

3075 

00.0 

1 1 

n 

4,04 

3,10 









Tube li: <1 = 

0,218 cm 









28 

40.34 

0.056 

14,08 

15,34 

4,28 

0,147 

K 

0.40 

0.009 

0.213 

411 

72.2 


012 

24,7 

yes 

no 

2.25 

37 

40.10 

15.42 

08.90 

114,32 

13.49 

0.133 

K 

0,03 

0.092 

0,304 

3005 

314,1 


5020 

70.9 

no 

yes 
• * 

2,89 

42 


15,40 

99,40 

114,80 

13,41 

0,133 

Ph 

0,00 

0,435 

0,10 

0.07 

0,414 

0,191 

3077 

277,3 

599,0 

2,10 

5030 

71.0 

•• 

3,03 

2,00 

43 


15.00 

90.75 

100.35 

14,07 

0.132 

Ph 

0,043 

0,420 

0.097 

0.000 

0.392 

0,174 

2848 

271,5 

010,0 

2.24 

4711 

08.0 

• i 

• i 

2,95 

1.93 

44 


10.33 

49,25 

05.58 

24,90 

0.118 

Ph 

0.715 

0.49 

0,095 

0.072 

0,427 

0.222 

1700 

153.0 

290.0 

1,93 

3250 

57,0 

•I 

ii 

3,28 

2,25 

45 


17.80 

41.15 

58,95 

30,15 

0.112 

Ph 

1,50 

1,03 

0.108 

0.072 

1,08 

0,400 

1580 

54.7 

120,0 

(2,31) 

3075 

55,5 

1 1 

ii 

7,29 

4,72 

40 


0.33 

15,15 

21,48 

29,40 

0,113 

Ph 

0,07 

0,048 

0,007 

0,055 

0,282 

0,224 

570 

70,2 

96,0 

1,20 

1110 

33.4 

yes 

no 

3,07 

2,97 

47 


0,50 

17,23 

23,79 

27,57 

0.115 

Ph 

0,578 

0,57 

0,005 

0,052 

0,250 

0,180 

038 

100,8 

128,0 

1,27 

1209 

34,8 

1 1 

n 

2,05 

2,01 

48 


4,82 

13,70 

18,52 

20,01 

0,117 

Ph 

0,333 

0,320 

0,001 

0,055 

0,127 

0,110 

497 

145,2 

167,5 

1,15 

925 

30,4 

ii 

ii 

1,53 

1,47 









Tube III: d = 

= 0,272 cm 









54 


85.5 

231.8 

317.3 

20.05 

0,110 

Ph 

2.050 

1,193 

0.194 
0.1 10 

2,50 

0,825 

5470 

127.0 

385.0 

3,03 

12850 

113 

no 

yes 

7,53 

4,39 

55 


80.5 

328.7 

415.2 

20,83 

0,123 

Ph 

1,500 

0.883 

0.180 

0,103 

1,754 

0,572 

7150 

237,0 

725,9 

3.00 

15790 

120 

• * 

ii 

5.51 

3.25 

50 


87,4 

350,1 

440,5 

19,57 

0,125 

Ph 

1,380 

0,815 

0,181 

0,080 

1,570 

0,441 

7090 

283,0 

1011 

3,58 

16700 

129 

•1 

ii 

5,10 

3,00 

57 


80.3 

374.3 

403,0 

19,20 

0.120 

Ph 

1,413 

0.800 

0,180 

0.088 

1,598 

0,443 

7990 

290,1 

1040 

3,61 

17300 

132 

n 

ii 

5,20 

2.94 


*) K = ca+he+ome+er Ph = photographed 
measurements 
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Figs. 1,2,3 



Figure 1.— Diagrammatic view of 
the test apparatus. 



Figure 2.- Apparatus 

for measur- 
ing the viscosity of 
gas mixtures. 



Figure 3.- Dynamic viscosity, 
kinematic viscosity 
and density for "propane" — 
oxygen mixture at 18 ± 1 °0. 
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Figs. 5,6,7 
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Figure 6.- Experimental values 

w fmax/ w fmln (solid 
curve) of the turbulent flame 
cone and comparison with the 

values (dashed curves) 

theoretically computed for the 
case of extremely fine-body 
turbulence. 


Figure 5.« Experimentally defined flame speeds Wf as 
function of the propane concentration. 




Figure 7.- Temperature variation in the laminar flame 
front at various hypothetical ignition 
temperatures and various f values. 
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Figs. 8,9,10,11 



Figure 8.- Expansion of a laminar free jet according to 
Schlichting. 



Figure 9.- Theoretically com- 
puted Bunsen cone 
in laminar flow. 



zzzzz luminous zone 0.12 tnjm thickness 

Figure 10.- Experimental Bunsen 
cone in laminar flow 
(test 47) compared with 
theoretically computed cone. 



Figure 

mixing 


■ Diagrammatic representation of a turbulent 
transport process according to the Prandtl 
length theorem. 
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Fige. 12,13,14 



Figure 12.- Turbulent exchange quantity e/uR in tube flow 
as function of the relative distance from tube 
center line r/R at several Re in dimensionless representa- 
tion. 



Figure 13.- Mixing length Z/R 
in tube flow as 
function of the relative 
distance from tube center line 
r/R for several Re in dimen- 
sionless representation. 



Figure 14.- Stream-line 
pattern of 
the turbulent marginal 
zone of a free jet 
(according to Tollmien). 
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Figs. 15,16,17 



nuclear and marginal zone of coarse-body turbulence 
of turbulent free jet at on an originally smooth 
various relative distances flame surface. 
x/R from the tube orifice. 
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Figure 17.- Flame cones computed theoretically for the 
extreme limiting case of fine-body 


turbulence 


